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Abstract 

We present experimental photoemission and inverse photoemission spectra of 
SrTi03_,5 representing electron doped d° systems. Photoemission spectra in 
presence of electron doping exhibit prominent features arising from electron 
correlation effects, while the inverse photoemssion spectra are dominated by 
spectral features explainable within single-particle approaches. We show that 
such a spectral evolution in chemically doped correlated systems is not com- 
patible with expectations based on Hubbard or any other similar model. We 
present a new theoretical approach taking into account the inhomogeneity of 
the real system which gives qualitatively different results compared to stan- 
dard homogeneous models and is in quantitative agreement with experiments. 
PACS numbers: 79.60.Bm, 79.60.Ht, 71.30.+h 
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Following the discovery of high T c oxides |T|, there has been a resurgence of interest 
in doped transition metal oxides. Recent investigations have revealed some very unusual 
physics in these systems [@-f§]. La x .Sri_a;Ti03 0j is a system containing Ti 3d electrons 
where the carrier concentration can be varied in a controlled fashioned. At one end, LaTiOs 
is a Mott-Hubbard insulator with a 3c? 1 configuration, while at the other end of the solid 
solution, SrTiOs is a band insulator with a 3d configuration. Thus, varying x tunes the 
electron configuration continuously between and 1. For large values of x, the transport and 
magnetic properties f2| as well as the experimentally observed spectral weights behave as a 
correlated doped Mott insulator. For small values of x (~ 0), the system represents small 
amount of electron doping in a 3g?° configuration, and thus, it should behave as a doped 
band insulator and is expected to show little correlation effects. This expectation is also 
consistent with transport and magnetic properties of such compounds which exhibit an 
essentially free-electron like behavior. 

Interestingly however, ultra-violet photoemission spectra (UPS) of La^Srx.^TiOs with 
small x (= 0.1 and 0.2) exhibit two features: a coherent low energy feature and a very 
prominent incoherent feature which is in fact found to be more intense than the coherent 
feature ||. The observation of a large fraction of incoherent to coherent spectral weight 
suggests the presence of very strong correlation effects, in contrast to transport and magnetic 
properties, as well as to the theoretical expectations of an uncorrelated behavior close to 
the d° configuration. Coherent features at low energies and incoherent features at around 
1.4 eV are also seen in systems with one d electron per unit cell [[|. The evolution of 
the redistribution of spectral weight as a function of interaction strength in the half filled 
system can be accounted by zero temperature large d calculations based on the Hubbard 
model || . The doping dependence of these features, however, has not been accounted within 
this framework. 

In order to clarify these important issues, we have reinvestigated the electronic structure 
of lightly doped d° systems both experimentally and theoretically. We supplement the al- 
ready existing UPS of La^Sri-^TiOa with those of SrTi03_,5, where the oxygen deficiency 
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(5) introduces electrons into the stoichiometric d° system. We have also carried out in- 
verse photoemission experiments in these systems. Since we expect that the strength of 
the incoherent features is more visible in the unoccupied part of the spectrum in systems 
which have a small number of d electrons, this technique is ideal for settling the importance 
of correlations here. The present results conclusively establish that the observed spectral 
function is inconsistent with the homogeneous Hubbard model which is normally used in 
the interpretation of these kinds of experiments. We present a new theoretical treatment 
of inhomogeneity (arising from disorder in the real system) within a Hubbard-like model to 
establish qualitatively new descriptions for the evolution of spectral functions in chemically 
substituted systems, in very good agreement with experimental results. Whereas substi- 
tuted compounds have been extensively investigated in recent times, the role of disorder 
has almost always been neglected; thus, our findings have important implications for a large 
number of compounds. 

Polycrystalline ingots of oxygen deficient SrTiC>3_5 with 5 = 0.0, 0.08 and 0.28 were 
prepared by melt quenching required amounts of Sr 2 Ti04, TiO and Ti02 in an inert gas 
arc furnace. Oxygen contents were determined thermogravimetrically. The sample sur- 
faces were cleaned by in-situ scraping using an alumina file at liquid nitrogen temperature 
to avoid surface degradation. UPS and bremsstrahlung isochromat (BI) spectra were ob- 
tained with 0.1 and 0.8 eV resolutions respectively, in a commercial VSW spectrometer [0. 
Self-consistent semi-relativistic ab-initio calculations with the linearized muffin-tin orbital 
method within the atomic sphere approximation (LMTO-ASA) for SrTi03_,5 with oxygen 
vacancies were carried out by considering a cubic supercell consisting of eight formula units 
of SrTiOs with empty spheres at oxygen vacancy sites and 216 /c-points in the irreducible 
part of the Brillouin zone. 

In Fig. 1 we show the experimental He i UPS of SrTi03_5 for 5 = 0.0, 0.08 and 0.28 near 
the Fermi energy (E^) region; the inset shows a typical spectrum over a wider energy range, 
being very similar for all the compositions with the intense doublet feature at about 7 eV 
arising from primarily oxygen 2p derived states. Expectedly, the wide band gap material, 
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SrTiOa shows no intensity at Ep. On introducing oxygen vacancies (5 > 0), it is clear that 
two spectral features emerge in the near region, one appearing at Ep and another at 
about 1.3 eV, in striking similarity to those of La x Sri_ x Ti03 [Q. While the feature at Ep 
is due to electron doping and the consequent movement of the Fermi energy Ep into the 
conduction band as confirmed by our band structure results on SrTi03_5, the calculations 
do not have any counterpart of the experimentally observed 1.3 eV feature (see later in the 
text). Furthermore, this higher binding energy feature in Fig. 1 (and also in La x Sri_ x Ti03 
0) is similar to the incoherent feature observed in LaTi03 |§ which is an end member of 
the solid solution La x Sri_ x Ti03. These observations along with the absence of this feature 
in single-particle calculations, would strongly suggest correlation effects being responsible 
for the origin of this peak in these lightly doped electron systems. 

In Fig. 2 we show BI spectra of SrTi03 and SrTi02.72 compared to the theoretically 
calculated BI spectrum of SrTiC>3 obtained from band structure calculations and calculated 
matrix elements |7j . The calculated spectrum provides a good description of the experimen- 
tally observed one for SrTiC>3 within the first 10 eV of E^; the mismatch at higher energies 
arises from the well-known limitation of linearized band structure methods. This agreement 
is consistent with the fact that SrTiC>3 with a d° configuration is primarily a band insulator 
with correlation effects playing no major role. However, it is interesting to note that the 
BI spectrum of SrTiC>2.72 is almost identical to that of SrTiC>3, except for a small shift 
0.5 eV) of the peaks and some overall broadening of the spectral features. The shift is 
consistent with electron doping of the system and thus, is explainable within a rigid shift 
of E F with respect to the one-electron results for SrTi0 3 with no perceptible signature of 
any correlation induced features in contrast to the occupied spectral features (Fig. 1). We, 
however, cannot rule out the presence of weak intensity features corresponding to the upper 
Hubbard band which may be considerably broadened if it coexists in the region of 4-9 eV 
with extended O and Sr states which have a much larger intensity. 

The above results of the presence of correlation induced feature in the occupied parts 
and the simultaneous absence of any such feature in the unoccupied parts, are very sur- 
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prising. It is well known that for few electrons in the band (such as Ce and U containing 
systems), correlation induced features appear most prominently in the unoccupied spectral 
function, whereas in Ni based systems with few holes in the band, the correlation effects 
are prominent only in the occupied parts. We demonstrate explicitly that the evolution 
of the photoemission and inverse photoemission spectra as a function of doping cannot be 
accounted by one band or a multiband homogeneous model by incorporating self-energy 
corrections to the calculated Ti d density of states (DOS) in SrTiC>3_5 (5 = 0.125) within a 
second order perturbative treatment of correlation [§]. The results for U = and 4.6 eV are 
shown in Fig. 3. It is clear from this figure that a correlation induced satellite appears out- 
side the one-electron (U = 0) bandwidth in the unoccupied parts with increasing U and no 
such distinct feature can be seen in the occupied parts. Similar results are also found within 
exact diagonalization studies of the Hubbard model || and in the solution of the Hubbard 
in infinite dimensions, while an exactly opposite trend is observed in the experimental results 
in Figs. 1 and 2. 

In order to understand this surprising result, we first note that the Hubbard model rep- 
resents a homogeneous system, whereas in real systems electron (or hole) doping is brought 
about by chemical substitution; this invariably introduces disorder or inhomogeneity at some 
length scale. In this paper we propose that it is essential to use a spatially inhomogeneous 
picture to model these systems. The basic idea is that in these systems the electrons experi- 
ence different local environments. The angle integrated photoemission spectrum, as a local 
probe, captures a suitable average of the Green's function over the local environments. In 
principle the degree of inhomogeneity can vary between two limits. If the system is disor- 
dered on a microscopic scale (termed microscopic inhomogeneity here) , but homogeneous on 
a macroscopic scale, i.e. the impurities do not form clusters, a model of uncorrelated disor- 
der is appropriate. If there is a tendency of the defects to form clusters, this may reinforce 
the natural tendency of strongly correlated systems to phase separate as pointed out in ref 
[IIJ and a model of macroscopic inhomogeneities is more suitable. Recent developments in 



the treatment of correlated electrons by using a mean field theory that becomes exact in the 



limit of infinite spatial coordination JTTJ] allows us to make these ideas more quantitative. 
We study the evolution of photoemission spectra in these two well defined limits. We expect 
that the physical situation is somewhere in between. We stress that in both limits we find 
the photoemission spectrum as a function of x to be very different from the predictions of 
the Hubbard model and very similar to what is observed experimentally. 

In the microscopically inhomogeneous situation we describe SrTi03_5 or La x Sri_ x Ti03 
using a Hubbard model with diagonal disorder. The Hamiltonian is given by: 

(ij) i » 

To get a bounded density of states and a well defined limit as the lattice coordination gets 
large, we take the hopping matrix elements to be defined on a Bethe lattice and scale 
them as -4^. The full bandwidth 2D = 4t, is estimated to be 2.5 eV from our band structure 
calculations. The £j are random variables which take the values €a with probability Pa and e# 
with probability Pb; Pa and Pb are given by the composition, as 25 and {1 — 25) in SrTi03_<5 
and x and (1 — x) in La x Sr 1 _ 2 .Ti03 respectively. Supercell calculations for SrTi0 2 .875 and 
SrTiC>2.75 showed that several bands are pulled down compared to SrTiC>3 due to the presence 
of vacancies by about 1.5 - 3.0 eV which provides an estimate for (e B — e A ). We use (e# — 64) 
~ 2.4 eV; however, the qualitative features of the calculations do not depend on this exact 
value. Ua = Ub = 4.6 eV was estimated from earlier spectroscopic analysis. The chemical 
potential \i is chosen so that the number of d electrons is equal to x in La^Sr^^TiOs or 25 
in SrTi03_5. 

The mean field theory which is exact in the d — > 00 limit determines the local Green's 
functions at the A and B sites as the local / electron Green's function of an Anderson 
impurity model 

H^J2^ata k( r + Uf+^f+f l 

k 

+ E V k (atU + f+a ko ) + e,f+f„ (1) 

k 

Here the / level position e/ is given by Ea — \i or Eb — n respectively. The hybridization 



function A(iu n ) 

has to be determined self consistently from the local Green's functions via the equation: 

A(iu n ) = t 2 G av (iuj n ) = t 2 (P A G A {iu n ) + Pb G B {ioJ n )) 

It is to be noticed that the A and B sites are connected by the same hybridization 
function which depends on the local physics at both sites. Thus we have a coupled system 
of equations. To solve the mean field equations we used the recent extension of the IPT 



method [12 1 to extract real frequency information. On the imaginary axis identical results 



are obtained with the exact diagonalization algorithm of Caffarel and Krauth ||13|| . 

We now turn to the case where the inhomogeneities occur on a macroscopic scale. In 
this scenario the density of electrons is inhomogeneous on a more macroscopic scale. A 
fraction fx is in an insulating phase with the local concentration of the Mott insulating 
phase riins = 1 and a fraction 1 — /i with the metallic concentration n met . The total density 
n = fin ins + (1 — fi)n met , and the measured Green's function is given by 

G = hG nins + (1 - h)G rimet (3) 

G nins and G„ met can again be calculated using the mean field theory now in the respective 
homogeneous phases. To explore this scenario we took n me t = .1 (which corresponds to a 
chemical potential \i = —.78) and varied the density n by varying the volume fraction fi 
in equation (|3]). The spectral function in these two scenarios for a density close to n = .2 
are shown in Fig. 4 a and b including a finite broadening to model the finite experimental 
resolutions. Fig. 4 a and b are in very good agreement with the experimental spectra (Figs. 
1 and 2); this is more striking so, when we keep in mind that we expect the real system to 
be somewhere in between the two limiting cases theoretically considered here. For example, 
the calculation correctly predicts the existence of an intense incoherent feature at about 
1.3 eV arising from the local Green's function for site A with a band-like coherent feature 



near Ep due to the local Green's function for site B. This is in perfect agreement with the 
experimentally observed features in the photoemission spectra of SrTiOs^ (Fig. 1) and 
LazSri^TiOs Moreover, it is clear from the calculated results that unoccupied parts of 
the spectral function will be entirely dominated by coherent band-like features arising from 
site B as observed in Fig. 2. Here we once again draw attention to the fact that these results 
are in sharp contrast to those expected on the basis of Hubbard and similar homogeneous 
models (see Fig. 3). While there have been past attempts to deal with substoichiometric 
correlated systems 114] , effects of disorder inevitably present in such systems have not been 
addressed so far. 

To summarize, in this paper we studied high-energy spectroscopic data of compounds 
with close to 3d configuration. We presented explicit calculations of the spectral function 
of correlated electrons at different dopings in an inhomogeneous picture. We considered 
two cases, in one the inhomogeneity is present only on a microscopic scale while in the 
other the system was assumed to be macroscopically inhomogeneous. Both models give 
similar angle integrated spectra which therefore cannot determine the length scale over 
which the system is inhomogeneous. The spectra of the inhomogeneous models, however 
are qualitatively different from that of a homogeneous system, conclusively establishing the 
need to go beyond the standard theories to describe spectral functions in a class of doped 
transition metal compounds. 
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FIGURES: 

Fig. 1 He i photoemission of SrTi0 3 _,5 for three values of 5. 

Fig. 2 Bremsstrahlung isochromat spectra of SrTi03 and SrTiC>2.72- 

Fig. 3 2nd order perturbation applied to the DOS obtained from LMTO supercell cal- 
culations for SrTi0 2 .875- 

Fig. 4 Spectral functions of the disordered Hubbard model for denisty n=.2 in the case 
of (a) microscopic and (b) macroscopic inhomogeneities. The separate contributions from 
the A and B sites to the total are also shown. The occupied parts of the spectral functions 
are expanded for clarity, as shown in the figures. 
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